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Abstract 



We describe and characterize a setup for subrecoil stimulated Raman spec- 

>: 

troscopy of cold cesium atoms. We study in particular the performances of 
^. ' a method designed to active control and stabilization of the magnetic fields 

o" 

across a cold-atom cloud inside a small vacuum cell. The performance of the 



setup is monitored by copropagative-beam stimulated Raman spectroscopy of 
a cold cesium sample. The root mean-square value of the residual magnetic 
field is 300 /j,G, with a compensation bandwidth of 500 Hz. The shape of 



the observed spectra is theoretically interpreted and compares very well to 



/\ • numerically generated spectra. 
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I. INTRODUCTION 

Stimulated Raman spectroscopy has become in recent years one of the most useful and 
powerful tools for laser manipulation of cold atoms. Raman stimulated transitions have 
been used, for example, in subrecoil laser cooling ||]; for the preparation of Bloch states in 
a stationary wave [Q]; to perform "sideband" cooling of trapped atoms ||; and in subrecoil- 
precision measurement of atomic velocity distributions in quantum chaos experiments [|J. 
This powerful technique is however plagued by its sensitivity to stray magnetic fields that 
should imperatively be compensated for. In the present paper, we describe a setup used 
in the above mentioned quantum chaos experiment for subrecoil-precision measurements of 
velocity distributions of cold atoms, with a special emphasis in a method for reducing stray 
magnetic fields to a sub-milligauss level by active compensation . 

All experiments reported below were done with cesium atoms, widely used in the present 
context, but most conclusions drawn in the present paper can be applied to other atomic 
species. Cesium recoil velocity v r (i.e. the velocity acquired by an atom emitting or absorbing 
a single photon) corresponding to the usual D2 line near 852 nm is 3.5 mm/s. The Doppler 
shift associated with a stimulated Raman transition at such a velocity is 8.4 kHz. 

Let us briefly review some relevant features of the Raman stimulated spectroscopy. In the 
present context, Raman stimulated transitions (Fig. [ft) are two-photon transitions connecting 
one hyperfine sublevel of the ground state to the other one by the absorption of a single 
photon in one Raman beam of frequency uj\ and wavenumber kx and by stimulated emission 
of a photon in the other beam (frequency uj 2 and wavenumber k 2 ) via a virtual excited level. 
The Raman process is resonant if the frequency difference of the two Raman beams and the 
ground-state hyperfine interval (ootf ~ 27rx 9.2 GHz for cesium) are equal, and the Raman 
detuning for zero magnetic field, low laser intensity and zero-velocity atoms is defined as 

5r = ujx-uj 2 ~u h f . (1) 

We suppose that the Raman beams are either parallel or antiparallel. In what follows we 
shall neglect \k\\ — \h%\ ~ \oJ\ — wi\jc compared to uii/c or a^/c, and take k = LOi/c ~ ujz/c. 



If the Raman beams are counterpropagating (i. e. k\ ~ — £2), the detuning seen by an atom 
of velocity v is 5r + 2k(v + v r ), and the atomic velocity sensitivity of the transition will be 
maximum and so is the momentum exchange between the Raman beams and the atom. This 
last property provides a way to perform subrecoil cooling |1|. On the other hand, if the two 
Raman beams have the same propagation direction (i. e. k\ ~ £2) the Raman detuning is 
almost insensitive to the atomic velocity As sensitivity to magnetic fields and to light-shifts 
is preserved, co-propagating Raman spectroscopy is very useful for the calibration of the 
apparatus. 

Since the life time of the F = 4 sublevel is very long (several thousand years in the absence 
of collisions), the width of the Raman stimulated transitions will be, in principle, limited by 
Heisenberg principle of uncertainty. The transitions linewidths are thus proportional to the 
inverse of the excitation duration. For instance, a Raman pulse duration of about 20 /xs, in 
the counter-propagating configuration, leads to a velocity resolution of v r . If the amplitude 
of the pulse is adjusted to produce a n pulse, it will transfer a whole velocity class of width of 
v r centered at the velocity Vq = 5n/2k + v r from one hyperfine sublevel to the other one. One 
can then reconstruct the atomic velocity distribution by measuring the transferred atomic 
population as a function of the Raman detuning 8r. Finally, in order to prevent resonant 
transitions to (and spontaneous emission from) the excited state, the virtual intermediate 
level must lie far enough from the closest excited level, what can be achieved by choosing 
an optical detuning A much larger than the linewidth of the closest excited level. Typically 
A = 200 GHz in our experiment, to be compared to the r = 2n x 5.3 MHz linewidth 
of the excited level. The spontaneous emission rate is thus reduced by about 9 orders of 
magnitude with respect to the resonance. This is true in the limit of a narrow laser line; we 
have however noticed that power diode lasers present a large spectral background about 40 
dB below the peak intensity that limits the reduction of the spontaneous emission rate. 

The optical part of our Raman setup has been described in detail elsewhere . Briefly, 
the Raman frequencies are generated by direct current modulation at 4.6 GHz of a diode 
laser. The current modulation generates two optical sidebands separated by 9.2 GHz that 



are used to perform injection-lock of two independent power diode lasers. We obtain in this 
way two 150 mW laser beams with the required frequency separation and with a beat-note 
width below 1 Hz. The Raman beams then pass through three acousto-optical modulators 
acting as optical switches that control the interaction time and allow shifting from the co- to 
the counter-propagating beam configuration. Finally, the resulting beams are transported 
by monomode, polarization maintaining, optical fibers to the region of interaction with the 
atomic cloud. 

II. ACTIVE COMPENSATION OF THE MAGNETIC FIELD 

The ground state of cesium splits into two hyperfine sublevels F = 3 and F = 4 separated 
by 9.2 GHz (the cesium-clock frequency). The linear Zeeman split for these states is of the 
form ZpBmp, where B is the magnetic field, nip the azimuthal quantum number and Zp 
the Zeeman coefficient (around 350 kHz/G in modulus) (Fig. |2|). An atom in the F = 4 
hyperfine sublevel undergoing a 6 mG magnetic field (2 % of the Earth's mean magnetic 
field) will experience a Zeeman broadening equal to the recoil velocity Doppler-shift. Thus, 
in order to perform atomic velocity manipulation with subrecoil resolution, stray magnetic 
fields should be reduced to the mG level or below. Usually, this is done by shielding the 
experimental zone with a double "mu-metal" sheet providing a high degree of isolation 
against magnetic fields. However, the shield reduces the optical access to the experiment 
and, as mu-metal is difficult to tool, it also tends to forbid further evolution of the setup. 
Furthermore, the magnetic gradient of a magneto-optical trap magnetizes the shield on the 
scale of a day can lead, leading to a slow magnetic field drift. It is thus interesting to 
perform an acceptable field reduction without any "physical" shield. In general the local 
magnetic field across the experiment is a sum of several components: DC earth field; DC 
stray fields, e. g. from an ion-pump permanent magnet; AC components from the line (50 
Hz or 60 Hz) and its harmonics; etc. This means that the magnetic field configuration is 
greatly dependent on the particular setup one wants to protect, and that a DC magnetic 



field compensation alone is insufficient. 

To overcome most of those problems, the method we describe here is to perform a real- 
time compensation of the magnetic field in the three directions with, for each direction, DC 
and AC compensations. In order to do so, we use the fact that the sum of the magnetic 
fields at the eight corners of a cube gives eight times the magnetic field at the center of this 
cube. This result is independent of the gradient of the field, and the first non-vanishing 
correction is of second order in the magnetic field |J 

In practice, we placed twenty-four magneto-resistive probes at the eight corners of the 
vacuum cell, three at each corner oriented according to the three orthogonal directions (see 
Fig. |3|). The probes are mounted on two U-shaped integrated circuits placed parallel to the 
top and bottom faces of the cell, preserving a large optical access to the cell. A maximum 
of algebraic operation is performed directly on the integrated circuits in order to reduce the 
number of external electric connections. The sum of the eight signals for a given direction 
generates three error signals: Each error signal is integrated and the result is used to drive a 
current supply (we simply use standard diode laser supplies in our setup) via a modulation 
input. Each supply drives a coil pair in Helmholtz configuration generating a magnetic 
field opposite to that measured by the probes. The rectangular compensating coils are 
typically 1.5 meter wide in order to insure the homogeneity of the compensated magnetic 
field across the vacuum cell (whose size is 10 cm). The magneto-resistive probes (HMC2003 
Honeywell) have, according to the manufacturer data, a resolution better than 40 fiG. In 
closed-loop operation, the residual magnetic field fluctuation measured by the probes is 150 
fiG, corresponding to a reduction of one order of magnitude, and the bandwidth is better 
than 500 Hz. We took care of aligning the probes axes with the axes of the compensation 
coils to better than five degrees, and we never saw any instabilities due to the coupling 
among the three directions with the three servo-loops acting simultaneously. With this setup, 
canning the DC magnetic field component along one direction (by scanning the respective DC 
current) does not affect the stability of the compensation along the other directions, even in 
presence of unavoidable stray couplings among directions, which are correctly compensated 



by the device. 

A good way of testing the effectiveness our setup is to perform co-propagating Raman 
spectroscopy. We have seen in Sec. | that such a configuration is almost insensitive to 
the atomic velocity, and it allows easy diagnostic of perturbing effects like residual stray 
magnetic fields and light-shifts. To do so, we programmed our experiment to execute the 
following time sequence: The magneto-optical (MOT) trap is loaded; the MOT gradient of 
magnetic field is shut down and a "Sisyphus" sequence is applied for further cooling. Twenty 
millisecond after cutting the MOT gradient, the magneto- resistive probes are reset. Reseting 
must be done before the probes can perform measurements once they have been submitted 
to intense magnetic fields (e. g. the MOT gradient). A special function is provided in the 
probe chip for this end. A few millisecond later, the servo loops are switched on. After 
twenty milliseconds, the optical beams are shut down, the atoms are transfered form the 
F = 4 to the F = 3 sublevel and the Raman pulse is applied. The Sisyphus sequence lasts 
for around 60 ms which is longer than in other cold atom experiments (typically 30 ms), but 
is much shorter than the requested delay (around 200 ms) when a mu-metal shield is used, 
necessary to allow for the shield magnetization relaxation. 

By applying a DC magnetic field during the experiment, one can perform spectroscopy of 
cesium ground-level Zeeman states. If the DC field has no particular orientation with respect 
to the polarizations of the Raman beams , the selection rules for the Raman transitions 
are Am = 0, ±1, ±2. As the Zeeman coefficient for F = 3 and F = 4 are almost equal and 
have opposite signs, Z 3 = 350 kHz/G and Z 4 = —351 kHz/G, the position (with respect to 
the line center) of the line connecting the sublevels (3, m) and (4, m') is 

Az/(m, m) (kHz) = [350(m + m) + m'} B(G) . (2) 

As we are dealing with weak magnetic fields, the second term in the brackets is in general 
negligible, and we can consider that the position of the lines depends only on m + ml. One 
easily deduces that the spectrum can display up to 15 lines, as shown in Fig. [|. This figure 
has been obtained by applying a DC field of 28 mG with the magnetic field compensation 



setup described above active. The Raman pulse duration is 500 /xs. In Fig. [| we compare 
a "zoom" of the line noted -6 in Fig. |] with the magnetic field active compensation on and 
off, clearly showing a noise reduction effect. When the compensation is off, only the DC 
component remains. 

The maximum measured full width at half maximum (FWHM) of the lines in Fig. ^ are 
around 1.5 kHz, which means that the resolution is independent of the magnetic sublevel 
and limited by the Fourier transform of the Raman pulse. For longer interaction times, 
we have observed a linewidth of 600 Hz which correspond to a velocity resolution better 
than i> r /10 and thus to a field residual amplitude of 300 fiG. To achieve this last value, we 
added a gradient-compensating coil pair in anti-Helmholtz configuration, that generates a 
10 niG/cm gradient to eliminate residual MOT-gradient magnetic fields. 

Reducing the DC field leads to a collapsed, Zeeman-compensated spectrum, with a 
FWHM of about 4 kHz (see Fig. |6|). In principle, one would expect the width of the 
collapsed spectrum to be just slightly greater than the widest line of the split spectrum. 
Explaining this "unexpected" broadening is the purpose of the next section. 

With the performances described above, we easily achieved a velocity-measurement res- 
olution of v r /2. As an example, Fig. ^ present a measurement of the velocity distribution of 
the cesium atoms. The FWHM is 82 kHz, corresponding to a temperature of 3.3 /iK which 
is typical of "Sisyphus-boosted" MOTs. This proves the ability of the setup to perform 
subrecoil spectroscopy. 

Concluding this section, we can say that this setup displays a sufficient resolution to 
be used in most Raman stimulated spectroscopy experiments. Furthermore, there is hardly 
any reduction of the optical access to the vacuum cell. Let us however mention that the 
setup presents a serious limitation in the everyday use. We observed that the HMC2003 
probes have a slow time-drift, with a time constant of roughly one hour. This implies that 
the DC component of the magnetic field compensation has to be readjusted regularly. We 
accounted for this effect by a software procedure restoring a fixed value of the magnetic 
field (as measured by the probes), for instance each hour. We also tested the possibility of 



replacing the magneto-resistive probes by small coils and use the induction signal to generate 
an AC error signal. This technique apparently presents no detectable drift. 

III. STUDY OF THE SPECTRA 

The relative intensity of the Raman transitions depends on the Raman amplitude con- 
necting the respective Zeeman sublevels, which is proportional to the product of the dipole 
matrix elements connecting each sublevel to the concerned intermediate excited level (see 
Appendix [A]) . The position of each Zeeman sublevel depends on the magnetic field and on 
the intensity of the Raman beams through the light-shift, which, for the far detuned Ra- 
man beams we are considering here, depends only on the dipole matrix element connecting 
the considered Zeeman sublevel to the excited level. We choose the magnetic field B to 
be parallel to the quantization axis, labeled z. The Raman beams are co-propagating with 
wavevector k (see Fig. [8]), and lie in the (x — z) plane, making an angle 9 with respect to 
the z-axis. The (orthogonal) polarization vectors associated with the two Raman beams lie 
in a plane orthogonal to k, and the polarization of the Raman beam starting from the level 
F = 3 makes an angle ip with the (x — z) plane (this last angle is irrelevant in the weak 
magnetic field limit, we shall thus always take if = in what follows). 

We consider three main causes of line broadening: (a) Magnetic field fluctuations and 
inhomogeneities, (b) fluctuations and inhomogeneities of the intensity of the Raman beams, 
and (c) the finite duration of the pulses. We shall consider each of these effects in detail. 

The fluctuations of the magnetic field displace the Zeeman sublevels, and thus displace 
the center of each Raman line (except the line connecting the Zeeman m = m' = sublevels, 
if allowed). Furthermore, we consider that the magnetic field in our setup results from the 
vectorial addition of a controlled, DC vector Bo and of a randomly fluctuating component 
AB. The resulting magnetic field thus varies both in modulus and in direction. The quan- 
tization axis is aligned with the DC component. The spatial inhomogeneity of the magnetic 
field also displaces the position of the Raman lines, thus contributing to the inhomogeneous 



8 



broadening. We do not know the exact distribution of the magnetic field intensity over the 
atomic cloud, but as the atoms have an almost random spatial distribution, this contribution 
can be statistically treated in the same way as a field fluctuation and we do not distinguish 
between them. The typical magnitude of this effect in the weak magnetic field regime is 1.5 
kHz. 

The dominant low-magnetic field broadening factor is the fluctuation of the light-shifts 
(as for the magnetic field effect, we do not distinguish the contributions of fluctuations and 
of inhomogeneities). These fluctuations are in part due to the fluctuations of the intensity 
of the Raman beams themselves. However, the light-shift affecting a given Zeeman sublevel 
(F, m) is independent of F and m if the direction of the magnetic field is parallel to the 
direction of propagation of the Raman beams (this result is also valid for zero magnetic field, 
as one can then always chose the quantization axis parallel to the beams). This implies that 
the light-shift fluctuations in this configuration displace the lines as a whole, and cannot 
contribute to the broadening of an individual line. This fact shall play an essential role in 
the following of this section. Thus, we do not expect the best compression of the spectrum to 
be achieved with zero DC field, because in this case the fluctuating component AB induces 
magnetic field random rotations with respect to k, leading to spectra that are broad and 
asymmetric. The thinnest lines are expected to be obtained with a DC field parallel to 
the Raman beams of an amplitude a few times greater than the mean amplitude of the 
fluctuating component: in such case, the fluctuation of the magnetic field induce only small 
misalignments with respect to propagation direction of the Raman beams. We verified these 
facts both experimentally and numerically (see below). Let us mention that a simple way 
to obtain the alignment of the DC field and the Raman beams is to try to minimize the 
amplitude of the lines corresponding to Am = ±1 (that is, -7,-5, ..,5,7) in the spectrum of 
Fig. |], as the Raman selection rules forbid such transitions if B and k are parallel. If the 
axes of the compensation coils are aligned with the Raman beams (as it is the case in our 
setup), one can reduce the DC component without loosing the correct field orientation. 

Lowering the intensity of the Raman beams produces thinner linewidths, but the over- 
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all signal also decreases, specially when performing counter-propagating (velocity-sensitive) 
Raman spectroscopy. The intensity used in our quantum chaos experiment [|J and in the 
curves shown in the present paper is 20 mW for each Raman beam of waist 4 mm. The 
typical value of the broadening associated with the light-shift fluctuations is 3.5 kHz. 

The last broadening effect is the finite duration of the Raman pulses. We worked with 
Raman pulses ranging from 1 to 3 ms for which the associated broadening (around 0.2 kHz) 
is negligible. 

Numerically calculated spectra have been generated by a Monte-Carlo technique. A 
particular value for the magnetic field is generated by adding a DC component Bo making 
an angle 9$ with k and a perturbation AB, whose amplitude is picked from a gaussian 
distribution of zero mean, root-mean-square amplitude AB , and of random orientation. 
Note that AB is formed by two contributions: the stray local fields (independent of B ) 
and the fluctuations of B itself, which we assume to be proportional to B Q (due e. g. to 
fluctuations of the current delivered by the supplies). 

In the same way, a value for the intensity of each Raman beam is picked from a gaussian 
distribution of mean intensity Iq and root-mean-square variation A/ - The position and 
the amplitude of each line is then calculated, and a particular spectrum generated using 
the "natural" linewidth given by the inverse of the pulse duration. Typically 1000 of such 
randomized spectra are averaged to produce the spectra shown in Figs. |] and |6|. 

As expected, the best collapsed spectrum is obtained for a non-zero value of the DC field. 
Fig. |9| displays the numerically calculated dependence of the FWHM on the DC magnetic 
field B Q , showing a well defined minimum around B = 0.5 mG. For stronger mean fields, 
the broadening is dominated by the fluctuations of Bq itself, and should be linear ion Bq. 
For values B close to zero, the broadening is due essentially to field rotations induced by 
the fluctuation component. The fluctuations of the field direction are roughly proportional 
to ABq/Bq and the related broadening is maximum if Bq = 0, but it cannot be infinite. We 
thus write the weak field term as A/(B + b), where A/b is the zero DC field broadening, 
depending only on ABq. The shape of the curve is then of the form 

10 



F ^=B^l +CB "- < 3) 

Fig. |^ shows that this forms fits rather well the numerical simulation. Unfortunately, we 
cannot directly compared this result to experimental data as we do not have an independent 
measurement of the mean magnetic field across the atomic cloud. 

Interestingly, one also finds that the minimum width does not correspond to 9 = (per- 
fect alignment of the magnetic field and the Raman beam wavevector) but to a small angle 
6q = 0.1 rad. For such an angle, light-shifts are slightly different for each Zeeman sublevel, 
and so are the Zeeman shifts. We attribute the minimum to a partial compensation of the 
two effects. However, this misaligned configuration tends to produce slightly asymmetric 
spectra. 

Concluding this section, we can say that we have provided a convincingly interpretation 
of the physical mechanisms governing the width of the Raman lines. The interpretation has 
been tested by comparing numerically generated spectra to the experimental results, which 
allowed us to quantify the performances of our setup. Moreover, this interpretation furnishes 
some interesting hints for the minimization of the width of the collapsed spectra. Another 
interesting hint is to polarize the atoms in one of the extreme (\m\ = F) Zeeman sub levels 
by optical pumping. This shall sensibly reduce the broadening due to fluctuations of the 
light-shift. Based on the measurements we have done in the split spectrum configuration, 
we expect in this case a collapsed linewidth of the order of 1 kHz. 

IV. CONCLUSION 

We presented in this paper an experimental setup for subrecoil stimulated Raman spec- 
troscopy. A key feature of the setup is the active compensation of magnetic fields over a 
small volume. Sensitive Raman stimulated transitions on cesium have been used to charac- 
terize the setup, with residual magnetic field fluctuations of 0.3 mG. The final limit of the 
method has been shown to be related to light-shift fluctuations rather than to the magnetic 
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field ones. The residual linewidth is 4 kHz, provided a small non-zero DC field is kept, al- 
lowing atomic velocity selection below vr/2 with the present method. Carefully constructed 
Raman 7r-pulses, would allow to reduce the light power, pushing the resolution on the atomic 
velocity selection to the ur/10 level, which would be comparable to the results obtained with 
double mu-metal shielding |2|]. 

Laboratoire de Physique des Lasers, Atonies et Molecules (PhLAM) is UMR 8523 du 
CNRS et de l'Universite des Sciences et Technologies de Lille. Centre d'Etudes et Recherches 
Lasers et Applications (CERLA) is supported by Ministere de la Recherche, Region Nord- 
Pas de Calais and Fonds Europeen de Developpement Economique des Regions (FEDER). 

APPENDIX A: LIGHT-SHIFTS AND TRANSITION RATE COEFFICIENTS 

The geometry and definitions used in this appendix are those of Fig. [8], explained in 
the preceding section. We choose the quantization axis parallel to the magnetic field and 
decompose the polarizations of the Raman beams in irreducible components (noted e_, eo, e+) 
with respect to that axis 0] 

— cos tp cos 9 + % sin tp 

£l = T2 

eio = cos tp sin 6 
cos tp cos 6 + i sin tp 
61 = 71 (A1) 

and the components of the polarization vector 62 of the second beam are obtained by chang- 
ing tp into 7r/2 — tp. 

The light-shift of the sub level (F, m) induced by the beam labeled i (i — 1,2) is propor- 
tional to 

L 4 (F, M) - (q Ei ) y Al (F,M;F',M') (A2) 

where q is the electron charge, E^ the electric field associated to the i Raman beam, p = 
{ — , 0, +}, r p is the irreducible p-component of the position vector, and (F 1 , m') correspond 
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to the Zeeman sublevels of the excited state (F' = 2, 3, 4 if F = 3 and F' = 3, 4, 5 if F = 4), 
and Aj(F, m, F', mf) is the optical detuning. For simplicity, we will consider that the optical 
detuning is high enough so that the above value can be considered to be independent of 
the sublevel (practically, this means that the optical detuning must be much larger than 
250 MHz). Applying the Wigner-Eckart theorem, one can reduce the above equation to the 
more useful form: 



Li(F,m) = (qE i ) 2 V2(2F+l)(2J+l)- 



\(QS 1/2 \\r\\QPj)\ 2 



E 



]Te jp v / 2F 7 TI 



1/2 J 1 



/ 



A 
F' 1 F 



\ 



m — p p —m 



(A3) 



F' F 7/2 

where the symbol into brackets is the reduced matrix element of r, the symbol into braces a 
"6j" symbol and the symbol into parenthesis a "3j" [§]. The total light-shift L(F, m) affecting 
the sublevel (F, m) is thus obtained by adding the contributions due to each Raman beam. 
Let us now consider the Raman transition amplitude (F 1; mi) — > (F 2 ,m 2 ), which is 
proportional to 



^^2^2(6S 1 / 2 F 1 m 1 \r p e lp \6P 3 /2F'm 1 - p) (6P 3 / 2 F'm 1 - p\r pl e 2p '\QS 1 /2F 2 mi - p - p) (A4) 

F' P p' 

that can be transformed into: 



|(6S' 1/2 ||r||6P J )| 2 (-l) i;,1+i72 ^- J + 1 /V2(2 J F 1 + 1)(2F 2 + 1)(2J+ 1 



^^ , 1 1 | _ 1/2 J 1 I I J 1/2 1 

J 1/2 1/2 \ f> \ F' Fi 7/2 \ \ F 2 F' 7/2 



p p 



I F> 1 fj W 



i mi — p p —mi 



\ 



F 2 1 F' 

, mi — p — p' p' — (wii — p) 



(A5) 
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FIG. 1. Raman stimulated transitions. 



m = -4 



m = _j, F=4(-351kHz/G) 

m = -2 



m = -1 



m = 

f m = 1 



m = 2 



m = 3 



OT = 4 



fflhf (« 9.2 GHz) 



OT = -1 



± m = 1 

. m = 



m = 1 



m = 3 



. m = -2 



m - " 3 F = 3 (350 kHz/G) 



FIG. 2. Hyperfine structure of the cesium atom ground state. 
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FIG. 3. Schematic representation of the probe geometry around the vacuum cell. The compen- 
sating coils are not shown. 
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FIG. 4. Experimental and numerically generated (see Sec. Il|) Zeeman spectra of Raman 
transitions in cesium. The lines of small intensity correspond to transitions Am = ±1, showing 
that the DC field is almost parallel to the Raman beams wave vector k. The numerical simulation 
was made with Bq = 28 mG, ABo = 0.5 mG, AIq/Iq = 0.025 and 9q = 0.2 rad. The asymmetry 
of the experimental spectrum is due to the fact that the Zeeman sublevels of the initial level are 
not equally populated. 
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FIG. 5. Detail of the line -6 of Fig. with the active compensation off and on. The non 
compensated curve presents a noise level close to 100%. The regular oscillations superimposed 
to the line itself are due to a stroboscopic effect between the network frequency (50 Hz) and the 
repetition rate of the measurements (4 Hz). 
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FIG. 6. Collapsed Raman spectrum at low magnetic field (FWHM 4.5 kHz). Comparison 
with the corresponding spectrum obtained by numerical simulation (see Sec. f||), obtained with 
B = 0.5 mG, AB = 0.3 mG, AJ /J0 = 0.027 and 6 = 0.01 rad. 
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FIG. 7. Velocity distribution of cooled cesium atoms. The FWHM is 82 kHz, corresponding to 
a temperature of 3.3 /iK, deduced from the gaussian fit. The resolution is of v r /2 (or, equivalently, 
4.2 kHz), showing the ability of the setup to perform subrecoil spectroscopy. 
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FIG. 8. Geometry of the Raman beams and the magnetic field. 
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FIG. 9. Theoretical FWHM of the collapsed spectrum as a function of the applied DC field. 
The fitting curve is 1.74/ (So + 0.21) + 3.24^0, with B in mG. 
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[6] Nevertheless, this is valid only for the exact geometric center of the cube. If the magnetic 
gradient is too high, the spatial extension of the atomic cloud can lead to an important 
broadening due to the inhomogeneity of the residual magnetic field. 

[7] We work with Raman beams of orthogonal polarizations, as, for beams far detuned from 
the optical resonance, the transition probability with parallel polarizations is zero. This 
is due to symmetry properties of the Clebsch-Gordan coefficients when summed of all 
hyperfine sublevels of the intermediate level 6P3/2. 

[8] Definitions and very useful formulas for calculation of these symbols can be found in the 
chapter 14 of Landau and Lifshitz, Mecanique Quantique, Theorie non relativiste, Mir, 
Moscou, 1966. Moreover, formulas (106,14), (108,10) and (106,10) are well adapted for 
the numerical computation of "3j", "6j" symbols and Clebsch- Gordon coefficients. 



